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I. INTRODUCTION. 

An important feature of hadronic collisions at high energies is the growth of the hard 
component of the interaction, induced by the increasing flux of partons. Apart from the 
single hard interaction, multi-parton hard interactions begin to play an increasing role. 
Hard events with multi-parton interactions have in fact been predicted long ago by 
several authors The simplest event of this kind, the double parton scattering, has 
been an object of the experimental search in all high energy hadron collider experiments 
since several years and while initially the results have been sparse and not very 
consistent, recently CDF has reported the observation of a large number of events with 
double parton collisions @. 

Multiple parton collisions represent a new observable feature of hadronic interac- 
tions. Non-perturbative inputs to the corresponding cross sections are the multi-parton 
distributions. These are new properties of the hadron structure which become acces- 
sible through the observation of multiple parton collisions. The multiple parton dis- 
tributions are related directly to the many-body parton correlations in the hadron [Q. 
While, as a general rule, the non-perturbative inputs are quantities to be measured 
and which cannot be computed in perturbation theory, in the interaction of two very 
virtual qq pairs, the whole process falls in the domain of perturbative QCD and, in 
that case, multiple parton distributions can be obtained through a direct calculation, 
within an approximation scheme. The aim of the present paper is to study precisely 
this example of high energy interaction and to work out explicitly the simplest case of 
multi-parton distribution. 

Note that studying events with several hard interactions one can introduce standard 
notions of inclusive and exclusive cross-sections. The first refer to the observation of 
at least one (or two,..) hard events, the second of exactly one (or two,...) hard events. 
The inclusive hard cross-sections are most easily accessable in the experimental study. 
They are also most easily calculated theoretically. In fact the inclusive cross sections 
are basically average quantities, since they include the multiplicity factor, which is 
proportional to the number of hard partonic interactions ||^, so no wonder that the 
inclusive cross-section may become larger than the total cross-section 0. The AGK 
rules tell that the inclusive cross-sections for single, double, triple etc hard collisions 
are given by the simplest relevant diagrams, which involve one, two, three etc hard 
interaction blobs. All other contributions from diagrams with more hard blobs cancel 
in the inclusive cross-sections. In the following, applying the simple theoretical formulas 
for the inclusive cross-sections, we shall use them as basic quantities, from which we 
shall extract the multiparton distributions. 

The paper is organized in two parts. For the sake of completeness the kinematics 
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of the single and double parton collision and the factorization of the non-perturbative 
parts of the processes are re-derived in the first part of the paper. In the second part, 
using the perturbative QCD approach, we compute the double parton distribution in 
two different limiting cases. The last paragraph is devoted to the conclusions. 

II. SINGLE AND DOUBLE PARTON INTERACTIONS 

A. Single hard scattering 

To fix our notations and to settle a common ground with the double parton colli- 
sion process, we discuss first the single hard scattering case which is described by the 
diagram shown in Fig.l. The discontinuity of the diagram through the hard blob rep- 
resents the inclusive cross section and, because of the AGK cancellation , it includes 
the factor accounting for the multiplicity of the hard partonic collisions in the hadronic 
interaction |^. In the following one considers only one sort of partons which, for the 
hard pomeron model, is represented by the gluons. The contribution of the diagram 
to the inclusive cross-section is related to the amplitude by 

a = (l/s)Im^ (1) 

where s = 2p+q_ is the cm energy squared. We neglect the masses of the projectile 
and target and take P- = q+ = p_l = q_L = 0. For massive projectile and target the 
standard substitution p ^ p — aq, q ^ q — ap with a = w? / s is assumed. The usual 
scaling variables of the colliding partons are introduces as 

fc+ = /_ = wq- (2) 

The partonic cm. energy squared is therefore written as Sp = = (k + l)"^ ~ 2fc+/_ = 
sxw. 

The standard treatment is based on two assumptions which are the basis for the 
factorization hypothesis 

* Parton virtualities (essentially the transverse momenta squared) are much smaller 
as compared with M^. 

** The dependence of the forward amplitude a on the virtualities of the external 
lines can be neglected. 

The longitudinal integrations in the diagram of Fig.|I] are easily done. Indeed ac- 
cording to our assumptions the amplitude a does not depend on A;_ nor on Z+. Only 
the upper part of the diagram depends therefore on /c_. One denotes it (with parton 
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legs) as Fi{p, k). Its integration over /c_ is standardly transformed into the integration 
over the "missing mass" si = {p — k)"^ of its right-hand discontinuity: 

/ ^F,ip,k) = — / ^lmF,ip,k) = F,ix,k^) 3 

Similarly for the target 

/f±F,(,.0 = ^ r ^ImF.(,,0 . (4) 

The standard parton distributions at the scale are obtained by integrating Fi 
and F2 over the transverse momenta kj_ and l±: 

p,ix,M') = J ^e{M'-e)F^{x,k) (5) 

The final integrations over and /_ can be transformed into integrations over x and 
w. The integration limits, [0, 1], are a consequence of the positivity constraints of si, S2 
and M^. The amplitude of Fig. 1 is therefore written as 

A^ = C '^-^p,{x,M')p2{w,M')a{M^) (6) 
JO xw 

while the single parton scattering inclusive cross-sections is expressed by the usual 
factorized formula: 

0"! = / dxdwpi{x, M'^)p2{w, M'^)ap{xw) (7) 
Jo 

The cross section can be expressed by using the unintegrated parton distributions in 
coordinate space. To this purpose one introduces the (non-forward) target parton 
distribution Fi{x, ki_i, k2±), shown in Fig.g, and its Fourier transform 

F,{x,n,r2) = J ^^^F,{xM,k2)e'''"'-''''' (8) 

(In the following part of this subsection one will not find 4-vectors, so we suppress the 
subindex _L for the /c's). It is convenient to introduce the partonic cm. and relative 
coordinates by 

Ti + r2 = 2R, ri — r2 = r (9) 

One has therefore 
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and when ki and k2 are equal 

= Fi{x,k,k) = J dFRd^rFi{x,R,r)e-''^ (11) 

which is precisely the unintegrated parton distribution which appeared in the single 
hard scattering expression above. According to (5) one has to integrate Fi over all k 
values, with the condition k^ « M^, in order to obtain the final partonic distribution 
at the scale M^. In the scaling approximation the p's do not depend on M^. One may 
therefore take the limit — » oo and thus integrate over all fc-values. In the scaling 
approximation one therefore finds 

pi(x, ^ cx)) = j (fRFi{x, R, 0) (12) 

When keeping into account the limits for k"^ the value r = is changed into ^ 1/M^. 
The standard parton distributions p{x, M^), entering in the factorization formula, can 
be therefore expressed via the density in coordinate space F as 

pix^M"^) = J d'^RF{x,R,r), = l/M^ (13) 

B. Double hard scattering 

The double hard scattering diagram is shown in Fig.|^. The corresponding inclusive 
cross section (including the multiplicity factor) is obtained, by means of the AGK rules 
0, by taking twice the imaginary part and dividing it by —s. 

The previous analysis can be repeated also in this case. There are five double 
longitudinal integrations. As independent variables one chooses ki,k2,h,l2 and q = 
ki — ks- The two hard scattering amplitudes are assumed to depend only on their total 
cm. energies squared = {ki + liY = 2ki^li_ and M| = (^2 + ^2)^ = 2^2+^2-, not on 
^1(2)- 5 ^1(2)+ nor on q±. The integrations on ^1(2)- are made as in the single scattering 
case. One calls the projectile blob (with all four partonic legs) Fi{p, ki, k2, ks, k^) 
{J2 ki = 0). Then the integrations over ^1(2)- can be transformed into integrations over 
the missing masses of the two hard scatterings Si = {p — ki)'^ and S2 = (p — ^2)^, the 
integrand being the discontinuity of F on the corresponding right-hand cuts. 

Since neither the hard scattering amplitudes nor F2 depend on q_ one may integrate 
Fi also over q_. This last integration can be transformed into an integration over the 
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missing mass in between the scatterings S3 = {p — q)^ along the right-hand cut, the 
integrand being the corresponding discontinuity. As a result one obtains 

^ dsi ds2 dss 



2p^ 



/ / 77-^/ 7r^'DisCs^BisCs^BisCs.,Fi{p, ki, k2, k3, h) 

Jsm 2m Js2o 2m Js-m 2m 



= -2tt i\ — Fi{xi, kij_, k3±\x2, k^,^, k2±) (14) 
\P+J 

The density defined on the right-hand side is the generalization of the single to the 
double unintegrated partonic density. The double parton distribution of the target is 
introduced in an analogous way. The integrations over ki(^2)+, h{2}- are transformed 
into integrations on the scaling variables Xi_2 and Wi^2 and the double hard scattering 
amplitude of Fig.| is expressed as 

i dxi dx2 dwi dw2 f -A- d'^kj d'^lj . ,3 

5'(fcl -k2-k3 + ki){27Tf6\ki +h-k3- k){2Tlf5\k2 + l2-ki- h) 

Fi{xi, ki, k3\x2, k4, k2)F2{x2, h, h\x2, U, l2)a'{Ml)a{Ml) (15) 

where Ml = sxiWi and M| = 3x2102- The double hard scattering cross section is 
therefore written as 



(T2 = / dxidx2dwidw2 / TT 



{2n)^ {2^f 

{2Tif5\k^ -k2-k3 + k4)i2n)^6\k, + h-k3- h){2TTf5\k2 + l2-k^- h) 

Fi{xi, ki, k3\x2, ki, k2)F2{wi, h, /3IW2, U, h)(^p{xiWi)(Tp{x2W2) (16) 

and both in Eq.(15) and in Eq.(16) the integrations are done with the limits kf, /c| < Mf 
and k1,kl< M|. 

The integrations over the transverse components of the transferred momenta are 
more conveniently performed by going to coordinates space. One introduces therefore 

Fi{xi,n,rs\x2,n,r2) = / U JTri^^^i^i^ ^i, k3\x2, k^, k2)e'^'''''~^-''''~^^^^^^ (17) 
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where in the integrand all four momenta are independent, in such a way that one is 
integrating also over the total momentum transferred to the projectile. One introduces 
the cm. and relative coordinates defined as: 

ri + r3 = 2Ri, n - = n^, r4 + r2 = 2i?2, ^2 - r4 = r24 (18) 

so that 

n = i?i + ri3/2, rs^Ri- ri3/2, r2 = i?2 + r24/2, r4 = i?2 - r24/2 
Similar coordinates, with primes, are introduced for the target: 

^1 + ^3 = 2i?i etc. 

Going to the coordinates space one finds, in Eq.(16), the following exponential 
exp I - iRi{ki - ks) - iri3{ki + k-i)/2 + iR2{k2 - k^) + ir2i{k2 + k^)/2 

-iR[ih - h) - ir[^{h + /3)/2 + iKik - k) + ir'2A{l2 + ^4)72} 

Also the three 5 functions in Eq.(16) can be represented as integrals over three addi- 
tional coordinates in transverse space: 

{2nf5^{ki -k2-k3 + k^) = J d2^e'ii(fei-fc2-fc3+fe4) 
{27rf 5\ki ^k-k^-h)^ j ci2^,^e*i(*^i+'i-*^3-'3) 

In the scaling limit the integrations over the eight momenta kj, Ij, j — 1, ...4 produce 
eight 5 functions: 

S^R + bi-Ri-ns), S\-R-b2 + R2 + r2A), 6\-R - h + Ri - ns), 
S\R + b2-R2 + r24), S\bi-R[-r[^), 6\-b2 + R'2 + 4^) , 
S^-b^ + Rl^-r'^^), S\b2-R'2 + r'2,) 
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The integrations over the transverse coordinates give then ri3 = = r24 = ^24 = 0, 
R[ = bi, R2 = b2, Ri = R + bi, R2 = R + b2. As a resuh the double hard cross-section 
is expressed as an integral over the three impact parameters, R, bi and 62 

(72 = / dxidx2dwidw2(yp{xiWi)ap{x2W2) 
Jo 

J d''Rd\d%2Fi{xi, R + 61, 0|x2, R + 62, 0)^2(^1, 61, 0\w2, 62, 0) (19) 

As in the case of the single parton scattering process, to keep into account that 
the integrations on the transverse momenta cannot be extended to infinity, one takes 
the relative transverse distances between the interacting partons to be of order of the 
dimensions of the two hard scattering blobs, namely is 1/Mi and I/M2. The resulting 
expression of the double parton scattering cross section is therefore written as 

(72 = y dxidx2dwidw2(Tp{xiWi)ap{x2W2) J d^ Rd'^bid^b2 

Fi{xi, R + bi,ri3\x2, R + b2,r2i)F2{wiM,rii\w2M,r24) (20) 

with rjg = I /Ml and rl^ = 1/M|. 

The expression is further simplified by introducing the cm. and relative parton 
coordinates in transverse space 

Bi = (61 + 62)72 + R, B2 = (&i + 62)72, 6 = 61 - 62 

The double parton distributions are therefore defined as 

Ti{xi,X2,b;Ml,M^) = J d'^BiFi{xi,R + bi,ri^\x2,R + b2,r2i) 

T2{wi,W2MMl,Ml) = j d^B2F2{wiM,riz\w2M,r2A) (21) 

and the double parton scattering cross section assumes a simpler form 

^2 = J Ti{xi, X2, b] Ml)ap{xiWi)ap{x2'W2)T 2{wi,W2, b; M^, M2)dxidwidx2dw2d% 

(22) 

with the geometrical interpretation shown in fig.^ Notice that the normalization con- 
tains also the multiplicity factor which counts the number of parton interactions in the 
hadronic collision (two in this case). 
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III. DOUBLE PARTON DISTRIBUTIONS IN THE BFKL REGIME 
A. Emission of a jet from a single BFKL chain 

In the framework of the hard pomeron model the hard colhsion process is associated 
to the production of a gluon with sufficiently large transverse momentum k > kmin and 
at a given rapidity (In this section only 2-dimensional transverse vectors will appear, 
so that we again omit the subindex ± and one take the metric Euclidean, k"^ > 0). 

The emission of a gluon with momentum k is described by modifying as follows the 
BFKL Green function in coordinates space G{r',r") [^]: 

G{r',r") J d^rG{r'r)V{k)G{r,r") (23) 

where r's are relative distances between the gluons and 

Vk{r) = ^ A e^'"- A (24) 

is the vertex operator describing the emission of the gluon, the arrows show the direction 
of the action of the A operators. The two BFKL Green functions in (23) are to be 
taken at appropriate energy (rapidity) ranges, corresponding to the rapidity distance 
of the emitted gluon from the projectile. This dependence is not written explicitly nor 
to overburden our formulas with arguments and/or indices. 

The first step is to generalize Eq.(23) for the non-forward BFKL chain. In momen- 
tum space the emission of a gluon with momentum k is described by the vertex 

Vihj2\l[J2) 

= ^i2nr5' (^(/i + h- l[ - Q - fc) {ill',' + lll[' - k\h - k?) (25) 
So the emission blob in momentum space Zk{qi, q'2|Q'i, Q'2) is given by 

z.i,u .^) ^ ^ / n 1111(2-)^^^ ili^^ + - - - ^) 

X (ill',' + lll[' - k\h - /2)')G(gi, g2|/i, /2)G(/i, I'M) 

X {2Txf5\q^ -q2-li + l2){2'Kf6\l'^ -I',- g', + g^) (26) 

In coordinates space one writes: 

(27r) V(gi - g2 - /i + /2)G(gi, gal/i, ^2) 

= / n^'^i^''2iG^(^i'^2ki,^2)e(-'''i''^+''^^"^+*"i-''^"^) (27) 
i=i 
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and 



i=l 

Note that due to the inclusion of the S function in the Fourier transforms (25) and 
(26), the Green functions in coordinate space are translationally invariant: 



G{ri + a,r2 + a\zi + a,Z2 + a) = G{ri, r2\zi, Z2) 



(29) 



The momenta squared in the emission vertex are substituted by differential operators 
applied to the Green function. The bracket in (24) is written as 

AiA^ + A2A; + k\Vi + V2)(v; + V^) 

where one uses the notation Ai = A(^i), A'^ = ■^(-^i) 6tc. The remaining 5 function 
in Eq.(24) is represented as an integral 

{2T:f5\{h + I2- l[ - l'2)/2 -k)^ j (i2^e-^^«'^+'^-''i-'^)/2-'=) 

in such a way that the integrations over the four transverse momenta in Eq.(24) produce 
four 5 functions in coordinate space: 

5\z^-z/2), 5\z2 + z/2), 5\z[-z/2), 5\z'^ + z/2) 

One obtains therefore 



,. 2 



i=l 



Akz 



X (AiA^ + A2A; + + V2)(V; + V^))G(zi, 4|r;, r',] 



^ ^/ z ^ 

Z\—Z^—-^,Z2—Z2- 



(30) 



and the differential operator which stands between the two Green functions is precisely 
the generalization of the emission operator \4 to the non-forward case. By using the 
notation: 

60;., 



r2 



G{ri,r2\zi,Z2)^e 



s Akz 



X (AiA^ + A2A; + A;2(Vi + V2)(V; + V2))G{z[, z'^\r[, r^) 



(31) 
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one may rewrite Eq.(28) in a more compact form 



ri + j(j2r-2 +iq[r[ - igj^j ) d'^ ^ 



xG(n,r2||,-|)v.(.)G(|,-||r;,rQ 
One defines Z in coordinates space in the foUowing way: 



(32) 



Zk{ri,r2\r[,r2] 



X {2'nf6{qr - q2 - q[ + q'2)Zk{qi, q2\q[, g;)e(^''^^-^'^^^-'''i'-i+^'^^^^) (33) 

By expressing the S function as an integral over the impact parameters B 

{27T)^6{qi -q2-q[ + q2) = J d^Be-'''^'^'-'^'~'''^+'^'^^ (34) 

and by using (32) and (34) in (33) one can do all integrations with the exception of 
those on B and z. The result is 

Zkiri,r2\r[,r'2) 

= J d'Bd'zG{n + B,r2 + b\^, -^)v,iz)G{^, -'-\r[ + B,r', + b) (35) 
and, as a function of the cm and relative coordinates, 



Zk{R,r\R', r') = J d^Bd'^zG{R + 5, r|0, z)Vk{z)G{{], z\R' + B, r') 
showing that Z^ is translationally invariant: 

Zk{R,r\R', r') = ZkiR + B, r\R' + B, r') = Zk{R - R', r, r') 



(36) 



(37) 



Eq. (37) is our final expression. It shows how one should change a BFKL Green 
function G{R,r\R' ,r') = G{R — R',r,r') to describe the emission of a gluon from the 
BFKL chain. 



B. Single parton distributions 



The simplest case to consider in the BFKL formalism is the inclusive emission 
of a gluon. The problem has been extensively studied in the literature [^]. Here 
we approach the problem from a somewhat novel point of view, introducing the single 
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parton densities in coordinates space in the BFKL framework. The quantities which are 
factorized in the BFKL formahsm are the 'unintegrated gluon densities' , depending 



exphcitly on the transverse momentum of the interacting gluon. The distributions 
that will be discussed hereafter are therefore different as compared to the distributions 
discussed in the first part of the paper, since they depend explicitly on the distance 
in transverse space between the initial and the final interacting partons. The latter 
distributions are therefore obtained only after integration on this transverse distance, 
with l//c as a lower limit. 

Rather than using the known expression for the inclusive cross-section in terms of 
forward BFKL Green functions, one starts form the amplitude, corresponding to the 
exchange of any number of pomerons between the projectile and target, obtained in the 
approximation of a large number of colors and assuming the coupling of pomerons 



to the colliding particles as purely perturbative |Tl[. The amplitude is written as 



A = 2is j (fR j (fr j dVpp(r)p,(r')(l - e-^^'^^W)^ (38) 

Here Pp{q){r) are the color densities of the projectile and of the target (with momenta 
p and q respectively). The function Gs{R,r,r') = G{R,r\0,r') is the BFKL Green 
function, depending on the Mandelstam invariant s. When retaining the s-wave only. 



because of the azimuthal symmetry, it takes the form [|T2 



Gs{R,r,r') = -— du ',,,,J d\J ^ (39) 

{2t:Y J [i/^ + l/AyJ Vriorso/ Vmo^2o/ 

where 

rio = R + r/2- Tq, r2o = R-r/2- Tq, r'^g = r'/2- Tq, r'^Q = -r'/2 - Tq 

and E{v) is the BFKL pomeron energy 

E{v) = (3//27r2)(ReV^(l/2 + tu) - V^(l)) (40) 

Eq.(38) allows us to write down easily the amplitudes for single, double etc. 
pomeron exchange processes. The single BFKL pomeron exchange amplitude is written 
as 

Ai = isgW^J d^R j d^r j d^r' pp{r)pq{r')Gs{R,r,r') (41) 

where the factor (7^ corresponds to the couplings of the pomeron to the color densities 
p. The inclusive cross-section for emission of a single gluon has therefore the form 
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d^R J d^r j d'^r'g'^pp{r)g^pg{r')d'^BdhGs,{R + B,r\0,z)Vk{z)Gs,{0,z\B,r') 

d^Rd'^R'd^z J rfV J d\'g'^pp{r)g^pg{r')Gs,{R,r\0,z)Vk{z)Gs,{0,z\R',r') (42) 

The dependence on energy is written explicitly in the Green functions. The emitted 
gluon rapidity (relative to the target) is y = logS2. Then si = s/s2 (the scale is 
not defined in the theory, as usual one takes it as 1 GeV). The integrated single hard 
scattering inclusive cross section can be written as 

0"! = y dyd^kl{k) = J d^k j dxdw6 ^ I{k) 

where one has introduced the scaling variables x = si/s and w = S2/S. By comparing 
the expression above with the factorization formula one obtains the expression for the 
partonic density of the projectile 

Fp(x, R, z)=g^ j d^rpp{r)Gs, (R, r, z) (43) 

with X = s/si and G{R, r, z) = G{R,r\0, z). 

The meaning of the variables is the following: R is the distance from the center of 
the projectile and z the distance between the initial and the final interacting partons. 
The size of z is of order of as a result of the factor ex^ikz in the hard interaction 
vertex, and the size of r is of the order of the dimensions of the projectile. The 
perturbative approach is justified when both r and z are small. The asymptotics of the 
Green function G for small relative distances and large energy was obtained in |TT| . 
The resulting expression is a function of the dimensionless variable ^ = R^ /rr'{» 1): 



_ / In^A 
Gs{R,r,r')r,r'«R^ I — Mn^exp — 



(44) 



where A is the BFKL intercept 



and 



A = ^ln2 (45) 

TT 
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42(1' 

a = ^C(3) (46) 

TT 

With z small and fixed, the partonic density (43) behaves as 

Fp{x,R,z) ^!/{^^^ -^logRexp ( ^^C^l (47) 

and it falls rather slowly (essentially as 1/-R^) as a function of the distance from the 
center of the projectile {R represents the distance from the projectile, whose average 
dimension tq = (r) is much smaller than R). The effect of the exponential is felt only 
at relatively large R ~ exp \Jl/x. A more quantitative estimate of the density can be 
made by neglecting the weak dependence on r in the logarithmic factors in (44) and 
by making the substitution r — rg: 

An analogous expression can be written for the target. 



C. Two pomerons coupling to a qq pair and triple pomeron 

The general structure of the amplitude describing the double partonic density in the 
BFKL regime is shown in Fig.|], where one assigns the partons undergoing the two hard 
collisions to two different BFKL pomeron ladders. The coupling of the two pomerons 
to the projectile is described by the upper blob B. One calls the energetic variable of 
the pomeron S2 = s/si, with s the overall cm. energy squared, for simplicity one takes 
it to be the same for both legs. The upper blob B is then integrated on si. 

The asymptotic behavior of B{si) is described by a single pomeron exchange, so 
that 

S(si) = B{y,) = 6(l/i)e^^\ yi = log s, (49) 

where A is the pomeron intercept and b{y) is some smooth function tending to a 
constant at high yi. Each of the pomeron legs has a similar behavior as a function of 
its energy variable S2'- 

P{S2) = P{y2) = P{y2)e^'^\ y2 = \ogS2 (50) 

with p{y2) a smooth function tending to a constant at high ?/2- The amplitude of Fig.§ 
is given by a integral over yi: 
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/ dy^B{y^)P^{Y - y,)P2{Y - y,) = e^^^ / dy^e~'^y%yMY - y^)p2{Y - y,) 
Jo Jo 

(51) 

where Y = yi + y2 = logs is the overall rapidity interval. As one may see, at very large 
Y values, the integration in (51) stays limited to the region where yi is smaller or of 
the order 1/A and thus it is independent on Y. One may therefore take for pi(2) the 
asymptotic values and rewrite (51) as 

^2AYpaSpas / ^-Ay^f^^^y^^^y^ (52) 
JO 

which corresponds precisely to the behavior of a two pomeron exchange process. 

The rapidities that dominate the integral (52) are of order 1/A. The actual value 
A is however not well defined, since it is given by the coupling, whose value is a 
parameter in the BFKL approach. Rather than trying to determine which are the 
rapidities relevant to B, we proceed by discussing two limiting cases. 

If is the virtuality of the interacting qq pairs (which has to be large to justify 
the perturbative approach) a limiting configuration is reached when is of the same 
order (or larger) than Si, namely, 

logg^ ~ (or >)1/A « Y (53) 

The blob B enters then in the DGLAP regime (finite scaling variable xi) where only 
the large logarithms log(Q^/A) are to be considered. In the fixed coupling approach 
of BFKL these logarithms are neglected altogether. Within the present scale invariant 
model one does not need therefore to sum any logarithms and one may state within 
the pure perturbative approach. The upper blob is then simply reduced to the qq 
pair, to which the pomerons are coupled directly. This is precisely the approximation 



which leads to the amplitude (36) obtained in |11] and which was used in the previous 
subsection. 

On the other hand, at lower virtualities, such as 

l«logQ^«l/A (54) 

the structure of the blob B becomes important. Another limiting configuration is 
reached when B enters in the BFKL regime, namely it is itself described by a pomeron, 
which has to eventually split into the two lower pomerons and the whole amplitude is 
described by a triple pomeron interaction. One should keep in mind that, in this last 
case, the three pomerons are not however in an equivalent regime. The two lower ones, 
as mentioned, may be taken at their asymptotical regime, Ay2 >> 1. The upper one. 
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on the contrary, is characterized by much smaller energies, such that At/i ~ 1, so that 
its exact form needs to be used. 

The two limiting possibilities which we consider are therefore the direct coupling of 
the pomerons to the constituent quarks and the triple pomeron interaction. We first 
discuss the simpler direct coupling case. 

D. Two pomerons coupled directly to the projectile(target) 

From the general expression of the amplitude (36) the contribution with two 
pomerons directly coupled both to the projectile and target is written as: 

g 

A2 = -is^- J (fR J (fr J d^r'pp{r)pg{r')Gl{R,r,r') (55) 

The integrated inclusive cross section is obtained by dividing the amplitude by s and 
by taking the imaginary part with a minus sign: 

02 = 9^ j d^R J d^r J d^r'pp{r)pg{r')Gl{R,r,r') (56) 

To obtain the double differential inclusive cross-section, corresponding to the emission 
of two hard gluons with momenta ki and k2 from the two pomerons, we have to 
substitute the pomeronic Green function with the emission functions Zk^ and Zk^- 
obtains therefore 

m„fc,) = i^^^a 
' dyid^ki dy2d'^k2 

g'J d^R J d\ J d\'pp{r)p,{r')Zk,{R,r,r')Zk,{R,r,r') (57) 
Explicitly, in terms of pomeron Green functions and emission vertices, one has 
I{ki,k2) = J d^Rj d\ J d^r'pp{r)pg{r') J d^ Rid^ R2dhidh2 

Gs,iR + Ri,r\0,zi)Vk,{zi)G,/,,{0, zi\Ri,r')G,,{R + R2,r\0,Z2)Vk,{z2)G,/,,{0,Z2\R2,r') 

(58) 

By comparing this expression with the standard form of the double hard scattering 
cross-section one may identify the double partonic density of the projectile 
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Fp{xi, Ri, Zi\x2, R2, Z2) = 9^^ j d?rpp{r)Gs^{Ri,r\{],zi)Gs^{R2,r\Q,Z2) (59) 

where logSi(2) = — logXi(2). A similar formula holds for the target. The density (59) 
can be further integrated over (l/2(i?i + R2) to obtain the double parton distribution 
depending only on the distance between the partons Vp{xi,X2, Ri — R2, zi, Z2) defined 
by (21). 

The expression (59) is a direct generalization of the single parton density (43) and it 
can be obviously generalized to the multiple partonic densities generated by the direct 
coupling of any number of pomerons to the projectile (target) 

/n 
(frpp{r)Y[GsXRhr\0,Zi), \ogXi = -\ogSi (60) 
i=i 

In this simple picture the multiple partonic distributions are factorized under the 
integral over r, which labels the different configurations of the qq pair. Interestingly 
this is precisely the correlation in transverse space which has been recently suggested 



T3[] to describe the anomalously small value of the effective cross section measured by 
the CDF experiment [0]. 

An estimate of the many-body parton distribution can be made, in this case, by 
using the asymptotics (42) and by substituting r in all slowly varying (logarithmic) 
factors by its (small) average value (r). One obtains 

Fp{xi, Ri, Zi\...Xn, Rn, Zn) ^ W Fp{Xi, Ri, Zi) (61) 

V / 1=1 

where Fp{x, R, z) is the single parton distribution defined by (41). Although the multi- 
parton distributions are just a product of single ones, they contain the factor (r")/ (r)" 
which is always greater than unity. The consequence is a positive correlations between 
partons, which however is independent of their position in transverse space. As an 
example the factor (r^)/ (r)^ is equal to A/tx for the Gaussian distribution p(r) and 3/2 
for the exponential one. 



E. Triple pomeron: generalities 

In Fig.p we show the amplitude in the case of triple pomeron interaction which, 
in the large number of colors limit, has been discussed in ref. |T^. The resulting 
expression is written as: 

g^N fvi , f d'^ri(fr2(fr3 

2 2 2 

r*-' /K*-" 
13' 12' 23 



Atp{si) = -isi—- I dy 
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^2; 2/1 - 2/)^'2(r2, ra; yi - y)r{^VlVl^!{ri, r^; y) (62) 

Here and in the future we use the notations ri2 — r-i — r2, R12 = {ri + r2)/2 etc. y's 
are the corresponding rapidity intervals. The pomeronic amphtudes ^ are defined as 
follows 

*(ri, ra, y) = j d^Kzd\[^P{r[^)Gs,{r[y^\n, r,) (63) 

with So = 6^ and similarly for the other two. To obtain the corresponding double 
partonic density, one has to substitute 1) the factor —isi by 2 and 2) the two pomeron 
amplitudes ^1 and ^2 by the appropriate pomeronic Green functions with free legs 
describing the gluons taking part in the hard interaction. The functions may have 
different rapidity intervals, corresponding to the different rapidities yi and 1/2 of the 
interacting gluons. The amplitude (62) is then converted into 



87r3 



dy j d'R[,d'r[,p{r[,) 



— 2 2% ^ Gs^is^iri, ral^i, ^3)^52/^0(^2, ^31^2, Z4)r\^VlVlG so{r[r'^\ri, rg) (64) 

^13^12^23 



where y = logso, Umin = mm{yi,y2}. 

Before moving further, one has to eliminate the differential operator acting on the 
last Green function, which makes the incoming pomeron unsymmctric with respect to 
the two outgoing ones in the triple vertex. To this purpose it is sufficient to note that 
the differential operator is proportional to the Casimir operator of the conformal group 

rts^lVl = I6C2 

Acting on a function 

/ ri3 y+^^" 

\r10r30J 

which appears in G{r[, r'^\ri, rs) it gives factor 16(i/^ + l/4)^. The action of the operator 
on G is therefore to transform it into a Green function which lacks the denominator in 
(39): 

a(i?, r, r')^-J duu's-^^^^ j d\o (^—) (^] ' (65) 

J J \rwr20J V ^10^20/ 
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Note that this function has the same asymptotics as G at high s, since it is determined 
by values of u close to zero. 

To obtain the double partonic density one has to finally take Z13 — Z24, — 0, and 
to fix the cm. coordinates of the initial gluons in both pomeronic legs R12 = Ri and 
R23 = R2- Note that the triple vertex coordinates i — 1,2,3 can be expressed via 
the cm. coordinates of the initial gluons 

n = R12 + R31 - -R23, r2 = R12 + R23 - Ri3, ra = R23 + R13 - R12 (66) 
so that 

ri2=2(i?i3-i?23), r23 = 2{Ru-Rl3): rsi = 2(i?23 - i?12) (67) 

The Jacobian of the transformation from ri, r2, to Ru, R23, R31 is 4. To obtain the 
density one has to drop the integrations over R12 and i?23- The double parton density 
in the triple pomeron interaction case is therefore 

g^N 1 rvmin f 2 f 

F(xi, Ri, Zi\x2, R2, Z2) = / dv d T p(r ) 

^ ' 1' II ^' ^' ^' UStt^ {Ri - R2)^ Jo ^ J '^^ ' 

I {R-R^\R-R2f^''''''^^'' ~ 

xGs,/so{R2, 2{Ri - R), Z2)Gso{R' - R, r', 2{R^ - R2)) (68) 

One can decouple the space integration involving the upper pomeron by shifting the 
integration variable R' ^ R' — R. The expression is then written as: 

g^N 1 

F{ X i,Ri,Zi\x2,R2,Z2) 



1287r3 - R2Y 
X / d''R!d\'p{r')Gs,{R:,r',2{R^-R2)) 

d^R 



X 



/ {R-R,r{R-R2r ^''^'' ~ R),zi)G,,/,,{R2, 2{R, - R), Z2) (69) 



F. Triple pomeron: calculation 

To evaluate the expression (69) one needs to make some simphfications: We use the 
asymptotic expressions for the two lower pomerons and we keep explicitly into account 
that the two relative distances zi^2 are much smaller as compared to the other distanced 
in (69), as a consequence of the large transverse momentum of the produced jets. A 
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further simplification is that the upper pomeron, although not in its asymptotic regime, 
is taken with zero total gluon momenta, where the Green function is substantially 
simpler. 

We first discuss the integration on it! (third hne in (69)): 

^'^^^^ I {R - rS{R - i^,)2 ^^iA°(^i' 2(i?2 - R), ^i)a./so(^2, 2{R, - R), Z2) (70) 

The asymptotics of the two Green functions at large 81^2/^0 and small 2:12 is easily 
worked out (see Appendix). The resulting expression is the same as in (42) with an 
additional factor 1/2 and a different definition of ^: 

G.,/.,(i?i,2(i?2 - i?),.i).,«i c ^;^H7^^-^er^ln6exp (^--^J (71) 

where 

\R — R\ — R2 1 1 ^ H~ -^1 — Ri I 

2\R-R2\z, ' y^-y^-y 

The same asymptotic expression holds for Gg^/so-i with 1^2. 

In the region 1/12 >> y >> 1 one may neglect the dependence on y in all smoothly 
varying factors and one may substitute yi 2 — > 1/1,2 everywhere, with the exception of 
the exponent. By shifting the integration variable R R — Ri — R2 one finds 

L(v) = r-2Ay Afai+ya) ^1^2 f 

^ ' n^a^ivmy/^ J R^\R + Ri\\R + R2\\R + 2Ri\\R+2R2\ 

In^iln^aexp (72) 



where 



R\R + 2Ri\ 
2zi\R + Ri\ 



, 6 = 6(1^2) (73) 



The integral (72) is dominated by the small R region, such that log(l/it!) ~ ^/yl^■ 
One then is allowed to put = in all factors which are finite in the R — hmit and 
Eq.(72) is simphfied to 



I^{y) = e-2AygA(2/i+j/2) 



■'1-2 



X 



47r^a^{yiy2r/' RlRl 

/ — - \n{R z^) \n{R Z2) exp 74 

J R^ \ ayi ay2 J 
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The integration over R now can be done explicitly (see Appendix) with the result 

L(v) = ^-2A?/ Afai+?/2) ( ^ ^1^2 



V «(l/l+Z/2) Z2J V «(?/l+l/2). 

At large yi and 2/2, with ?/i = 0(2/2)5 the second term is sub-dominant and it may be 
dropped. The expression is therefore reduced to 

h{y) = e-^^e^(-^-)4 ( ^^V' 4%exp ( J^M^] (76) 

The next step to evaluate (69) is to integrate the forward Green function, appearing 
in the first line, over its rapidity variable with the factor arising from Ii. 

roo 

h = / dye~''''Gs„Q,,=o{r', 2(i?i - R^)) (77) 
Jo 

Here it has been explicitly indicated that the Green function has to be taken at total 



momentum of the gluons equal to zero. One then obtains |12 



^,(0, r', r) = — / rfz/s-^(^)(r/r')''" (78) 

where -E'(z^) is given by (40). By putting this expression into (77) and by integrating 
over y one is left with an integral over u 

dv{r/r'f (79) 



27r2 7 ' ' ' 2A + E{u) 

where r = 2(i?i — R2). The distance r' is of the oder of the projectile dimension and 
it is small, so that the ratio r/r' for fixed -Ri — -R2 is a large number. The integral 
(79) can be therefore evaluated by taking the residues of the integrand at the zeros of 
2A + E{v) in the upper half-plane. 

The zeros are located at the points v = ixk, Xi < X2,< .... The first three points 
are (see jlSl ) 

xi = 0.3169, X2 = 1.3718, x^ = 2.3867 (80) 



One obtains 
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/2 = ^rrXcfe(r/r')-'"'= (81) 



where 



c^^ = V'(l/2 - Xk) - ^'{1/2 + Xk) (82) 
At large values oir/r' the nearest pole contributes, so that one finds 

h ~ -^c^rr'ir/r')-'^''^ (83) 

OLgN 

The final average, namely the integration on r' with the color density of the projec- 
tile, gives therefore (r^-^^), to be compared with (r^) which is the result obtained for 
the double parton density when coupling the two pomerons directly to the qq pair. 

Collecting all the factors the double parton density corresponding to the triple 
pomeron picture is finally written as 

F{xi,Ri,Zi\x2,R2,Z2) 



Z1Z2 / \yi'{Zi/z2) \ 

RlRl\R,-R2V+^^-^''^\ a{y,+y2)) 

(84) 

Analogously to the case of direct coupling of the pomerons to the qq pair, the double 
parton density contains the factor 1/R\R\. However in the triple pomeron case the 
factor — i?2|^^^^^ induces an additional strong positive correlation in transverse 

space between the two partons. 



IV. CONCLUSIONS 



In the present paper wc have evaluated the double parton distributions, in the case 
of interactions between very virtual qq pairs in the BFKL regime. After factorizing the 
hard interactions, the double parton distributions depend on the structure of the blob 
attached to the projectile (or to the target). Two hmiting possibihties have therefore 
been considered: a) the whole interaction is represented by the exchange of two BFKL 
pomerons attached directly to the qq pair, h) The process is described by triple pomeron 
interaction, where the blob is itself represented by a BFKL pomeron. The two-body 
correlation in the double distribution are different in the two cases. In the first case 
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the only correlation in the double parton distribution is induced by the configuration 
taken by the qq pair in transverse space, namely at a given qq configuration the resulting 
many-body parton distribution is just a Poissonian. The second case is, on the contrary, 
characterized by a strong correlation in transverse space, as the common source of the 

two partons is a BFKL ladder rather than the qq pair. The two possibilities are linked, 
through the AGK rules, to the diffractive events and to the events with rapidity gaps. 
The relative importance of the two contributions can be therefore inferred by comparing 
the rates of the two processes. 

APPENDIX 

In this appendix we shall give some technical details concerning the asymptotics 
(79) and calculation of the integral in (72). 

First the asymptotics of Gs{R,r,r') at small r'. From (37) one concludes that at 
s ^ oo small values of u contribute. So one is lead to study the leading behavior of 
the integral over tq at small u when r' 0. In this limit the second factor in (37) is 
singular at Tq = 0. So small values of Tq give the dominant contribution and one can 
take the first factor out of the integral over vq at ro = 0. The integral over ro then 
simplifies to 



(To keep the integral convergent one should take Imi/ > 0). This integral is trivially 
calculated by going over to the momentum space: if 




(85) 




(86) 



then 




(87) 



Function f{q) is trivially found by the inverse Fourier transformation: 



f{q) = 7r2^+2-5-^-^-^ 




(88) 



Putting this into (83) and doing the integration we find 



/ = (i7r/2)r' 



(89) 
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Inserting this expression into the integral over v together with all accompanying factors 
and taking its asymptotics at s — > oo one obtains the desired asymptotics (79). 

As to the integral which appears in (72), in the variable (3 — \nR\t reduces to an 
integral of the Gaussian type 



which can be done in a straightforward manner, with a result (73). 
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FIG. 2. Soft blob in the single scattering term. 



FIG. 3. Double scattering term. 




FIG. 4. Graphical representation of Eq.22 
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FIG. 5. Double parton density in the case of direct coupling of the BFKL pomerons to 
the constinuent quarks. 




FIG. 6. Double parton density in the triple pomeron case. 



